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Abstract 

The differential cross-sections for elastic and inelastic neutrino-deuteron scat- 
tering are calculated analytically using nucleon-nucleon effective field theory. 
For elastic scattering, the deuteron axial form factor and the deuteron strange 
magnetic moment form factor are computed to next-to-leading order, includ- 
ing two-body currents. For inelastic scattering, two neutral current processes 
vd unp, Vd — > Vnp and one charged current process Vd e^nn are com- 
puted to next-to-leading order. These depend on an isovector axial two-body 
matrix element whose value is yet to be fixed by experiment. Potential model 
calculations by Kubodera et al. and Ying et al. are reproduced for differ- 
ent values of the two-body matrix element. This implies that the differences 
between the two potential model calculations lie in their treatment of short 
distance physics. The charged current to neutral current vd cross-section ratio 
is confirmed to be insensitive to short distance physics, and the same ratio is 
obtained by potential models and the calculation presented here, within 5%, 
for incident incident neutrino energies up to 20 MeV. The two-body matrix 
element could be fixed using the parity violating process ltd enp. 
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I. INTRODUCTION 



The Sudbury Neutrino Observatory (SNO) employs inelastic neutrino-deuteron scatter- 
ing to study the solar neutrino flux, and seeks the solution to the solar neutrino problem. 
The cross-sections for these processes have never been measured, leading to a reliance on 
theoretical calculations in the analysis of their data. To date, fairly sophisticated potential 
model calculations [0,0, employing different approaches to the inclusion of meson exchange 
currents have produced results which agree to within 5-10% of one another up to neutrino 
energies of 160 MeV. Still, how well constrained are these calculations? Given the impor- 
tance of SNO to the understanding of neutrino physics, these cross-sections merit critical 
study. 

Ellis and Bahcall first calculated the cross-section for the reaction 

Ue + d^e^+p + p (1) 

in 1968 10, to aid in an early attempt to build a deuterium solar neutrino detector. Later 
work by other authors focused on using the reaction 

Ue + d—^Ue + n + p (2) 

as a probe for neutral currents and a test of various electroweak models of the time, in 
particular the Glashow-Salam- Weinberg model [^]. 

The SNO proposal for a large scale heavy-water solar neutrino detector led to a 
resurgence of interest in both processes, and to the need for precise theoretical calculations 
of the cross-sections up through supernovae neutrino energies of order 20 MeV, and beyond. 
While there have been a number of calculations of progressive complexity |[^-|5|,|9|-|TI| , there 
are two definitive calculations of these cross-sections which include modern nucleon-nucleon 
potentials and meson-exchange current effects, though the approaches are quite different 
IHJ^. The differences between these calculations are quite small, and it would appear that 
the results have converged. 

However, there are no definitive experimental tests of these calculations, yet SNO must 
rely on them in order to extract information on the solar neutrino flux. As such, it is 
important that the physics behind these calculations is understood thoroughly, along with 
any underlying theoretical and systematic uncertainties. 

Our approach to the problem of neutrino-deuteron scattering is to take advantage of 
recent developments in the use of low-energy effective field theory (EFT) |]T2|-^^. The 



power counting scheme of Kaplan, Savage, and Wise [21| allows for a concise and systematic 



analysis, order by order in perturbation theory, and has been used to study dynamical 



processes involving the deuteron, including electromagnetic form-factors and moments [23 



Compton scattering [|2^j25|] , np d'~f |2^j2^, and pp — > de~^i>e, where electromagnetic 



effects are important [^. It is straightforward to analyze neutrino-deuteron scattering in 
this scheme. First, we will look at elastic u-d scattering because of its sensitivity to the 
strangeness property of the deuteron. Then, we will proceed to two processes of interest to 
SNO, the neutral current (NC) reaction u + d —>■ u + n + p and the charged current (CC) 
reaction i? + d^e^ + n + n. The normal CC reaction, u^ + d e~ + p + p, is not considered 
here because of the additional complication of electromagnetic effects in the two nucleon 
final state. 
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II. THE LAGRANGIAN 



A. Effective Field Theory 

The lagrangian for an effective field theory involving nucleons and mesons can be de- 
scribed via 

C = Co + Ci + C2 + --- , (3) 

where contains operators involving n nucleons. Neglecting, for the moment, the weak- 
interaction couplings 

A = ^TrD^.ED'^St + ^ATrm,(S + St) + • • • , (4) 
o 4 

where E is the conventional unitary representation of the pion fields in SU(2) 



E.exp^^j , n = ^-r_;^^j . (5) 

/ = 132 MeV is the pion decay constant, and the trace of the quark mass matrix rUq is 
related to the pion mass through = \{mu + rrid)- The single nucleon lagrangian is 



+ 



(6) 



where the a Pauli matrix acts on nucleon spin space and ^ = exp(in//). The nucleon axial 
coupling is qa = 1.26. 

The two nucleon lagrangian needed for next-to-leading order (NLO) calculations is 



ci 



'2 



{N^PiNy{N^(D^Pi - 2& • + Pi^^)N) + h.c. 



{ 'So) 
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8 



{N^PiN)\N^(D^Pi - 2& • Pi^ + Pi^'^)N) + h.c] , (7) 



where Pi and Pj are spin-isospin projectors for the ^5*1 channel and the ^5*0 channel, respec- 
tively, with definition and normalization 
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Pi=^-a2T2n , TrPjp,- = 1(5,,- , (8) 

where the r matrices act on isospin indices. The strong coupling constants Cq, C2 and D2 
have renormalization scale (/x) dependence. The details of the fitting procedure for these 
parameters, along with their values in each channel, can be found in the Appendix. 
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B. Weak Interactions 



The effective lagrangians for charged and neutral current weak interactions, in terms of 
neutrino, nucleon, and meson fields, are given by 



--^/^J'.^h.c. , (9) 



= -^/zJ', , (10) 
where the is the leptonic current and is the hadronic current. For u-d and V-d scattering, 

= 1^7^^(1-75)6 , Fz^vn^-^,)u . (11) 
The hadronic currents can be decomposed into vector and axial-vector contributions 

j-^v--A; = ivi-Al)-tiv^-Al) , 

J^^-2sm'ewV^ + {l-2sm^ew)V;^-A'^-Al , (12) 

where the superscripts represent isovector components (with S representing isoscalar terms) 
and, later, the currents will be labeled by the number of nucleons involved. 

In a next-to-leading order (NLO) calculation, the electron mass mg contributions to 
the matrix elements are counted as higher order, along with contributions to the current 
proportional to the momentum transfer since q^l^ = up to NLO. Weak-couplings to pion 
fields are also higher order and are neglected here. The non-relativistic one-body isoscalar 
currents are given by 
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Ar = -'-AsN^^^^^N 
° 2 2Mn 



Af'^~AsN^a,N . (13) 



Similarly, the isovector currents are given by 



'k 



2 

2Mn 



N 



Af^gAN^^a.N . (14) 
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There are vector currents from magnetic moment terms, with k^^^ = ^{np + and k*^^^ = 
^{Kp—Kn) being the isoscalar and isovector nucleon magnetic moments in nuclear magnetons, 
with 

Kp = 2.79285 , K„ = -1.91304 . (15) 

The isoscalar weak currents are also sensitive to strangeness matrix elements between nucleon 
states. As measures the strange quark contribution to the proton spin 

2Sf,As = {p |s7^75s| p) (16) 

with the value determined by Savage and Walden 

As = -0.17 ±0.17 (17) 

from an analysis of lepton scattering data [^,0 including SU(3) symmetry breaking effects 
| 47| . is the covariant spin vector, fig is the strange magnetic moment of the proton 



{p |s7^s| p) = Up{G%{q ) + G'^iq 



2Mn ' ^ ' 



G^(0) = , iis^Gim . (18) 



The value of [ig has a large uncertainty. The SAMPLE experiment |^9| has measured 

G'i^(-O.lGeV^) = 0.23 ±0.37 ±0.15 ±0.19 n.m. , (19) 

However, theoretical predictions for /i^ range from —0.8 n.m. to 0.18 n.m., as summarized 
in ref . fS^ . 

Finally, there are two two-body axial currents relevant to scattering at NLO: 
Af = -2iei,kL2AN''m)\N''P,N) , 

Af = L,^A {N'^PkNf (N^PaN) + h.c. . (20) 

In addition, there is a two-body vector current that contributes to the NLO strange magnetic 
form factor 

Vk^^'^ =2zmN^PiNy(Vi + Vi){N'^PkN) + h.c. . (21) 



III. NEUTRINO-DEUTERON AND ANTINEUTRINO-DEUTERON ELASTIC 

SCATTERING 

A. Deuteron Form Factors 

ud ud and Vd Vd elastic scattering processes are sensitive to the strange quark 
properties of the deuteron. This can be easily seen by writing the neutral current in terms 
of quark degrees of freedom. 
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117^(1 -75)M-c?7^(l-75)rf-s7;.(l-75)s -2sm'^wJT' , (22) 
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then taking the matrix element between deuteron states. The fact that the deuteron is 
isoscalar means that the matrix element depends only on the electromagnetic and strangeness 
properties of the deuteron, 



( d J' d) 



d 



-2 sin' ^wJ"J" 



1 



d) + ( d 



1 



d 



d 



d)-( d 



A1 



d 



(23) 



Since the deuteron's electromagnetic properties are well known, ud and ud elastic scattering 
processes are, in theory, ideal for studying the strangeness of the deuteron and might also 
provide valuable information on strangeness in the nucleon. 

To begin, we must define the NC vector and axial form factors of the deuteron, then 
calculate the strange form factors in EFT up to NLO. Using \p,i) to represent a deuteron 
with momentum p and polarization state i, the vector form factors of the neutral weak 
current can be defined in the way similar to that of the electromagnetic form factors 



P ,3 



1 



2M, 
1 

+ 2M| 



2Mi 

[F|(q2)%(p + p'), + F!M'){5-^^ - 



(24) 



(25) 



where Md is the deuteron mass and q = p' — P- The neutral weak charge, magnetic and 
quadrupole form factors F^, F^j and Fq are linear combinations of the electromagnetic 
charge, magnetic and quadrupole form factors Fc, Fm and Fq and the strange charge, 
magnetic and quadrupole form factors F^, F^j and Fq. 



fS = 


-2 sin' 


^wFc- 


-F6 

2 ^ 


Fm = 


-2 sin' 


'^wFm — 




FS = 


-2 sin' 


^wFq- 





(26) 



The electromagnetic form factors of the deuteron are well known experimentally, and have 
been computed in using effective field theory. The form factor Fq = at q' = 0, because 
there are no net strange quarks in the deuteron. This means that F^, , together with Fq do 
not contribute to the cross-section until next-next-to-leading order (NNLO) and will not be 
considered in this paper. 
The axial form factor 



P ,J 
p',J 



A! 



iEjkiFl^G^) 



2^( P + P )fc 
2M, 



(27) 
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arises directly from strange matrix elements only, as has been noted from eq.(p3D already. 

Fl= FX . (28) 

Note again that terms proportional to have been dropped in eq.(p7D since we treat the 
neutrino as massless. 

The effective field theory result for the strange form factors and F|j up to NLO can 
be written as a series expansion in powers of Q, where Q is the small expansion parameter 
in the momentum and chiral expansion 

ps ^ pm ^ ps(i) ^ . . . _ ^29) 

The superscripts denote the power of Q in the expansion. The LO strange axial form factor 
depends on the strange axial moment of the nucleon As (defined in eq.([T6|)) and the LO 



electric form factor F(j 



(0) 



-AsFS^iq') , (30) 
where 

frtf)^^ta„-.(M) . (31) 



7 = y/M]sfB, and B = 2.2245 MeV is the deuteron binding energy. 

The NLO strange axial form factor depends on the NLO electric form factor F(j \ the 
electric quadrupole moment of the first non- vanishing order Fq~^\ and a two-body operator 
with coefficient L2^a from eq. (|20|) . 



Ff\e) = -As(F«(q^) + ^F^^^(q^)) + L,,M^^^-^ . (32) 



12M]^^ (q)) + ^2,A(/i)^^ 



-^2,a(/^) itself depends on fi but the form factor is /i-independent. The expressions for Fj^^ 
and Fn can be found in 



The strange magnetic form factor is computed as 

fg"(q') - 4^.(4"(q') + + mp) "'"'^''^ - , (33) 

where Hs is the strange magnetic moment of the nucleon defined in eq.(18) and L| is the 
coefficient of the two-body counter term defined in eq. (pT]) . 

B. Cross-Section 

The i'{y)d elastic scattering cross-section up to NLO is 
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da 

dVt u(V)d 



27r2 



cos^ - (4 sin* •§ 



w 



+2 sm^ - 



(34) 



where u is the incident neutrino (antineutrino) energy, and 6 is the scattering angle between 
the incident and outgoing lepton directions. The negative sign of the magnetic and axial 
interference term corresponds to ud scattering and the positive sign to Vd. Note that in 
addition to the dependence which is not considered in we also disagree on the sign 
of the interference term. In our expression, a^d < <^ud when F^^ = 0. This can be confirmed, 
later, when we consider the elastic limit of the inelastic scattering process. 
In reference 



the motivation was to measure As through the ratio 



R. 



■el 



daiyd ud) — da{ud ud) 
daiyd ud) + da{ud — ^ ud) 



(35) 



At LO, this ratio depends on As and Hs while at NLO it also depends on L2,a and Lj. 
Given that the values of As and fig should become known more precisely from single-nucleon 
experiments, precision measurements of the ratio Rei could tell us more about the intrinsic 
strangeness of the deuteron. 



IV. NEUTRINO-DEUTERON AND ANTINEUTRINO DEUTERON NEUTRAL 

CURRENT INELASTIC SCATTERING 

A. Neutral Current Structure Factors 



For neutrino scattering from a hadronic target, the differential cross-section can be writ- 
ten in terms of leptonic and hadronic tensors l^i, and W^i, as 



d^a G?.|k'| 



duj'dn 327r2|k| 



(36) 



where k(k') represents the initial (final) lepton three-momentum. The leptonic tensor is 
given by 



r = 8{k>'k'" + k^k"' - k ■ k'g"" + ie>"'f"'kpkl) 



(37) 



The hadronic tensor can be defined as the imaginary part of the forward matrix element 
of the time-ordered product of two weak current operators. It can be parameterized by six 
different structure functions 



W, 



J dW'^^T (c/(P) \jf{x)J^{0)\ d{P) 



— Im 

TT 

P P P'^qf^ 



Ml 



Ml 



(38) 
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where the momentum transfer = k^ — k'^ (the difference between the incident and outgoing 
lepton four-momenta), and is the deuteron four-momentum. W^, and Wq do not 
contribute to the differential cross-section because qj.^^ = 0. 
For the reaction 



u + d^u + n + p, 

the differential cross-section in the lab frame (deuteron rest frame) simplifies to 



duj' dVt 



27r2 



2Wi sin^ - + W2 cos' 



2 M. 



W3 sin^ 



(39) 



(40) 



where 9 is the angle between k and k , cu' is the final lepton energy, and we have used the 
relation 



q^ = — 4cc;u;'sin' 
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(41) 



For ud — > Vnp scattering, the last terms on the right hand sides of eq. (pTf ) and (|40|) change 
sign. 

The phase space boundaries for this reaction through the angular bound 



Max 



-1,1 



4:Mm(u - B) - 



< cos 6* < 1 



and the bound on outgoing lepton energy 



Q<uj' <uj- 2{Mn - ^Ml - 72) 
where v = uj — uj' is the energy transfer (go)- 



(42) 



(43) 



B. Leading order 

The leading order contributions to the hadronic tensor W^^y are the diagrams shown in 
Fig. |l]. The non-zero components of W^y are given by 



J0)2 



Wl^<^ = 5.,{2[cT + &f)F, + \{Cr 



{If 



)F2 + AdfF, 



C 



(1)2 



F2 + ^d!^~F. + ^C 



(1)2 



(44) 



with no summation implied over indexes of Wl(^ . The coefficients arise from the specific 
forms of the single-nucleon vector and axial currents, and are given by 



C 



(0) 



(0) 



— sin^ 

-Us 

2 



w 



C 



(1) 



1 



(1) 



9a ■ 



(45) 
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+ 




+ 




+ 



FIG. 1. Leading order contributions to the v-d scattering hadronic tensor. The gray blobs arise 
from no insertion plus any numbers of insertions of four nucleon operators Cq or Cg as 
appropriate. The crossed circles denote operators that create or annihilate two nucleons with the 
quantum numbers of the deuteron. The solid lines are nucleons. The solid circles correspond to 
insertions of leading weak current operators. 

The functions Fa are the individual contributions at LO from the diagrams of Fig. |l|, 
with F3 and F4 both associated with the last diagram in the first row: 



Fi = Re 
F2 = Re 



2MnJ P 



71 [M2,z/2 - qV] 

^tanh- fM^^ 
7rz/|q| XMj^u 



2 7 



|q| 



TT |q| 

F, = F,{^S,^ 'So) 



2(7 -ip) 



(46) 



The magnitude of the relative momentum between the final state proton and neutron is p, 
with 



P = ]l Mnu - 7^ - — + ze 



(47) 



The 'rescattering' amplitudes appearing in F3 and F4 are the leading order NN scattering 
amplitudes in each channel, given by 



A'y.'\p) 



1 + 



(48) 



(/i + ip) 



and similarly for the ^5*0 channel. After matching onto the effective range expansion as 
performed in the Appendix, 



A^lt'\p) 



An 



1 



M, 



N 



(49) 
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and 



A 



An 1 



M 



Nl + ip 



(50) 



where a*-^"^"^ is the ^5*0 channel scattering length. 

From these expressions for the components of the hadronic tensor, we find that the 
leading contributions to the structure factors in eq.(^OD are 



W. 
W. 



LO 



' * '11 



LO 



LO 



LO 



LO 
00 



LO 



(51) 



1. Elastic Limit 

An important test of the inelastic calculation is that we can reproduce the elastic scat- 
tering results in the limit 

p — ^7 + e . (52) 

It is straightforward to see that, in this limit, 

Fi = , 
^2 = , 

F4 = . (53) 

This, combined with eq.(|y,|^,0), reproduces the LO contribution to eg. (|3^ . Further it 
can be shown that our NLO result also recovers the elastic limit. 



2. Threshold Behavior 

Another useful test of the inelastic calculation is to study the threshold behavior and 
compare that with the results expected from the effective range expansion. In the effective 
range expansion, it is well known that the dominant contribution to the threshold hadronic 
matrix element is the ^Si ^Sq transition through the isovector axial coupling. The 
^5*1 — > ^5*1 transition is suppressed because amongst the NC spin-isospin operators 1, r", cr* 
and r"(T*: i) the isovector operators don't contribute (the transition is isoscalar); and ii) the 
matrix elements of the isoscalar operators vanish in the zero recoil limit {d and np states 
are orthogonal in the zero recoil limit). 

Our results reproduce these features. In the threshold limit ,p — > 0, 

VVqq — u , 

W,f°= (M;^) (!-««■=») 7)^ . (54) 
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FIG. 2. Graphs from insertions of the operators with coefficients D2 "^^^ or D2 '^"^ (denoted by 
the gray squares) that contribute to the v-d scattering hadronic tensor at NLO. The gray blobs and 
other features are as defined in Fig. 

Another consequence of the ^Si — > ^Si suppression is that our results will not be sensitive 
to isoscalar parameters As, fig and L2^a- This has the added advantage of reducing the 
number of free parameters, as will be seen later. 

Using eqs.(|y,^,0), the differential cross-section with respect to the relative kinetic 
energy Er{= /Mn) between final state nucleons is 



— — = / rf COSt^ ZTT-— — — 

dEr J-1 ^ ' du'dVl 
^ 2Glg\{u - B fMMP 



(1 _ a('^^o) _ 



This reproduces the threshold behavior of the effective range expansion result of [^]. Our 
result is also consistent, at both LO and NLO, with the analytic expression of the np d'j 



amplitude given in . 



C. Next to Leading Order 
At NLO we will decompose the hadronic tensor into five components 

w^,'^'' = w;^^ + w;:,- + w;:j + w^: + w^:f- , (56) 

each corresponding to an insertion of a NLO C2 operator, single pion exchange or higher 
order weak couplings. Each of these five components will be considered separately. 

1. D2 Contributions 

The D2 contributions to the hadronic tensor arise from the diagrams shown in Fig. |^. 
They can be expressed as simple replacements of the LO results in eqs. (|i^ , |46| ), 

(57) 
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FIG. 3. Graphs from insertions of the two-body current operators with coefficient Li^a (denoted 
by the circles) and L2,A fffroy circles) that contribute to the v-d scattering hadronic tensor at NLO. 
The gray blobs and other features are as defined in Fig. |||. 

where 

and similarly for the ^5*0 channel. 



(58) 



2. La Contributions 

Diagrams with one insertion of Li ^ or L2,a contribute to the hadronic tensor at NLO 
are shown in Fig. ^. These contribute only to the spatial piece of the hadronic tensor, and 
we find 



Svr^ |q| I \2(7 — ip) 



J(59) 



with all other = 0. 

Note that the term involving Li^a behaves like 

which is not /i-independent. The (3 function of Li^a can be obtained by requiring the 
A^A^ —>■ NNuV, ^Si to ^5*0 transition amplitude to be /i-independent, 




^(i5o)^(35i) 



(61) 



3. C2 Contributions 

Fig. ^ shows the diagrams in which C2 contributes to the hadronic tensor at NLO. Here 
we find contributions to both spatial and time-like components of the hadronic tensor 
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<30 OO 



FIG. 4. Graphs from insertions of the operator with coefficient C2 or C2 " (denoted by the 
solid squares) that contribute to the v-d scattering hadronic tensor at NLO. The gray blobs and 
other features are as defined in Fig. |i|. 



|q| 



00 



TT^ |q| 



Im 



ALf'Hp) tan 



-1 



.2(7 - w), 



(62) 



and 



W 



C2 



Im 



v4L'^^(p) tan 



-1 



|('5o) 

0,C2 

q| 



2(7 - ip)^ 



37r2|q| 



AL'f°)(p) tan-i 



|q| 



.2(7 -^p). 



(63) 



where 



^0,C2 



(P) 



(64) 



and similarly for the ^Sq channel. All other W^^ = 0. An analysis using eq.(|6TD shows that 
the sum of C2 and Li^^ contributions to the hadronic tensor are /i-independent. 



4. Single Potential Pion Exchange Contributions 

The potential-pion exchange diagrams at NLO do not yield closed-form analytical solu- 
tions. However, a reasonable set of approximations could be introduced to make analytic 
results possible. The diagrams contributing at NLO are shown in Fig. ^. For all diagrams 
q dependence is neglected in the pion propagator. In the region < 20 MeV that we are 
interested in, |q| < 40 MeV and the error due to this approximation is estimated to be 
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FIG. 5. Graphs from potential pion exchange that contribute to the v-d scattering hadronic 
tensor at NLO. The dashed lines are pions. The gray blobs and other features are as defined in 
Fig. |. 

(q^) /m^ < 10%. Further, for the diagrams in the second row of Fig. we angle-averaged 
the pion propagator; 



(p — k)2 + ml (p2 + k^) + ml 

where p and k denote the nucleon momenta at the tcNN vertices. This approximation could 
lead to a 30% error at this order, but numerically the error is found to be of order 10%. 
Since the terms neglected here are formally of order NLO, we have not performed a complete 
calculation to NLO. But the error due to the approximations made accumulates to no more 
than 20% of the NLO contribution from potential pions, and is numerically an NNLO effect. 
The structure factors can be decomposed in a manner similar to the LO result, with 

WS = 2{df + df )G, + liGf - G^f)G, + iGf 'G, + |C« . (66) 

The components are suppressed by factors of (q^) /ml compared to the diagonal 

ones, so we neglect them (if they contribute at all). All other = 0. The functions Gi^2 
are simple modifications of the LO functions Fi 2, given by 

G, = FJ{p) , G2 = F,f{p) , (67) 
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where 



m^r + 27 p 



- tan 



P 



(68) 



G3 and G4 can be further decomposed into 3 parts: 



Lt4 — -|- -|- LT4 



(69) 



with 



=i^3(A^T^(p)-4:r(p)) , 



1(^51) 



G. 



(a) 



(70) 



where 



_ 9\ ( m^M^AL'f )(p) 



2/2 



47r 



+ 



(71) 



(6) _ n^b) 



tan-f^^Utan-^l^l+2^^ 



27 



27 



+ mi tan ^ 



2(7 + <^m2+p2jj 



(72) 



Finally, 



7r|q|/2 i 

1 ^ 



^vr^ |q| 



tan 



|qi 



^2 (7 - ip) J + 27 



+ li + ip] + /i(7, p) + h{-ip, ij) 



tan 



.2(7-^p), 



(73) 



Gr = Gi^'i^S, ^ 'So) , 
where 



h{j,p) 



1 m^ + -f 
H — ^ in 



167r^(m^ + 7) i ^ — ip 7^ +p' 
q2 + 4(72+p2 + |q|p) 



m^r + 27 
m^r + 7 - ip 



|q| 



-q2 + 4(72 +p2 |q| 7) 



(74) 
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FIG. 6. Graphs contributing to the v-d scattering form factor at NLO. The dark solid 
circles correspond to insertions of the leading weak current operators, while the gray solid circles 
correspond to insertions of the next-to-leading weak magnetic moment operators. The gray blobs 
and other features are as defined in Fig. |i|. 

5. Contributions 

The diagrams with one LO single-nucleon weak current and one NLO weak magnetic 
moment couphng, as shown in Fig. ^, contribute to the antisymmetric part of the NLO 
hadronic tensor 

= -^W^j^^eo.,,^ , (75) 
with all other = 0. It is straightforward to show that 

-fc^^^C^^F.-lc^^CffF, , (76) 

where 

c£) = -2sin2^?^^KW-^/i5 , Cff = (l-2sin2^H')«:« . (77) 

Recall that at LO, vanished. As such, this NLO contribution is the leading contribution 
to the difference between neutrino and antineutrino scattering from the deuteron. 



V. ANTINEUTRINO-DEUTERON CHARGED CURRENT INELASTIC 

SCATTERING 



The differential cross-section for the process 

z7 + (i— i>e^ + n + n 



(78) 



takes the same form as eq. (^). There are, in principle, corrections due to the finite positron 
mass me, but their impact on the total cross-section is negligible, even at threshold (they 
could, however, affect the angular distribution). We must also modify the phase space 
bounds to include the effects of the positron's mass, as well as the neutron-proton mass 
splitting 6m = nir, 



nip. The modified bounds are 



Max 



-1. 



4Mjv(z/ - B - 6m) - u;^ - cj'^ + m^ 



2uJuj'^ — ml 



< cosO < 1 



me<uj' <uj- 2{Mn - ^ M% - Mn{B + 6m)) . 



(79) 
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For the most part, however, the primary difference between the neutral current and 
charged current cases is the fact that the charged current processes are purely isovector. 
As a result, the charged current structure factors can be obtained from the neutral current 
structure factors with the simple substitutions: 





= , 


Lyy — 




= , 


_ 

*-"A — 


L2,A 


= , 
= 


— 

Li^A — 



IK 



ud 



\Vud\ 



-9a 



V2 

z/ ^ z/ - (5m , (80) 



where we use IV^^I = 0.975 for this CKM matrix element. 



VI. NUMERICAL RESULTS 

(Pa 

We have presented analytic expressions of the differential cross-sections , for vd — >■ 

duj dil 

i>np, ud — > unp and ud — > e~^nn processes in the previous sections. These expressions have 



the correct elastic scattering limits and their threshold behaviors are consistent with p8 
We can now consider the numerical results for the total cross-sections of these processes. 
In calculating cross-sections for the different channels, we will include the isospin splitting 
(charge dependence) of the strong interaction in the ^Sq channel, but it can be neglected in 
the weak interaction currents as a higher-order effect. 

In the NLO calculations, four parameters (As, /i^, -Z^i,a, and I>2,a), which are not well 
constrained, contribute with different significance. To demonstrate their numerical effect we 
express the NC cross-sections at = 10 MeV as 



a{u{u)d u{u)np) = 0.999 ± 0.026 + O.OlSLi,^ n 
+ 10-5As(±0.5 ± 1.2/i, + 6.3As - 4.6^2^ ^ ^ 



where Li^a and L2,a are in units of fm^. It is clear, immediately, that As and fig contribute 
less than 1% to the total cross-section. Furthermore, since dimensional analysis suggests 
Li A and L2,a are of order 

\Li,a\ ^ \L2,a\ ~ ~ 5 fm^ (82) 

M 

at yU = m^r, we see that Li a could contribute at the 10% level but that L2^a would contribute 
far less than 1% to the total cross-section. This agrees with the ^5*1 ^Si suppression 



mentioned in section [IVB 2| . In this computation, Li ^ is kept as a free parameter. As = 
—0.17, /is = 0, and L2,a = 0. The error from this choice is estimated to be less than 1%. 

The LO results are parameter free, and are shown in Fig. ^ together with potential model 
results. With respect to the two most modern potential model calculations used, Kubodera 



18 




7.5 10 12.5 15 17.5 20 

E^{MeV) 




7.5 10 12.5 15 17.5 20 

Ep{MeV) 




8 10 12 14 16 18 20 

E,{MeV) 



FIG. 7. Cross-sections for vd unp (top graph), vd unp (middle graph) and vd — > e'^nn 
(bottom graph) shown as functions of incident energy. The solid curves correspond to the 
average of YHH and KN results while the dashed curves correspond to the LO results in EFT. Note 
that at LO, there is no difference in EFT results for ud — > vnp and ud — > Dnp cross-sections. 

et a/.'s results (KN are persistently larger than Ying et al.'s results (YHH |l|]) by 5-10% 
in all channels. Thus, we average the KN and YHH calculations in each channel, before 



^We note that while these calculations are presented in Kubodera and Nozawa of ref. [Q], they are 
actually the unpublished work of Kohyama and Kubodera, also listed in ref. [|| . 
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^.(p) MeV 


ud - 

a 


vnp 

b 


vd — > 

a 


unp 

b 


ud - 

a 


h 


4 


0.0288 


0.000299 


0.0283 


0.000299 






5 


0.0888 


0.000965 


0.0869 


0.000965 


0.0264 


0.000268 


6 


0.188 


0.00212 


0.183 


0.00212 


0.111 


0.00120 


7 


0.330 


0.00381 


0.319 


0.00381 


0.262 


0.00299 


8 


0.516 


0.00609 


0.496 


0.00609 


0.484 


0.00576 


9 


0.747 


0.00899 


0.714 


0.00899 


0.780 


0.00960 


10 


1.02 


0.0125 


0.973 


0.0125 


1.15 


0.0146 


11 


1.35 


0.0167 


1.27 


0.0167 


1.59 


0.0207 


12 


1.72 


0.0216 


1.61 


0.0216 


2.11 


0.0281 


13 


2.14 


0.0271 


2.00 


0.0271 


2.70 


0.0368 


14 


2.61 


0.0334 


2.42 


0.0334 


3.35 


0.0467 


15 


3.12 


0.0404 


2.88 


0.0404 


4.08 


0.0580 


16 


3.69 


0.0480 


3.38 


0.0480 


4.88 


0.0706 


17 


4.30 


0.0564 


3.92 


0.0564 


5.74 


0.0846 


18 


4.97 


0.0656 


4.50 


0.0656 


6.68 


0.0999 


19 


5.68 


0.0754 


5.12 


0.0754 


7.68 


0.117 


20 


6.45 


0.0860 


5.77 


0.0860 


8.74 


0.135 


TABLE I. 
in the form a 


NLO ud and ud inelastic 
= (a + hLi^A) X lO-^^cina 


scattering results in EFT with the cross-sections written 
with Li^A in units of fm^. 



comparing them to our LO results. First, consider the NC channels. At LO there is no 
difference in our calculation between ud unp and Vd — Vnp cross-sections. We see then 
that our LO result lies within 10% of the (averaged) potential model results. A similar level 
of agreement is seen in the comparison of our LO calculation of I^fi —>■ e^nn to (the averaged) 
potential model calculations. 

Our cross-sections at NLO will have one free parameter, Li^a- To facilitate study of the 
effects of this two-body contribution, we tabulate the NLO cross-sections as functions of 
Li^A in Table |. It is not surprising that we can reproduce the total cross-section results 
of KN and YHH in all three channels to a high degree of accuracy by choosing Li ^ to be 
6.3 fm^ and 1.0 fm'^ respectively (see Fig. |^). This would imply that the 5-10% systematic 
errors in these two potential model calculations are due solely to different assumptions made 
about the short distance physics. It also suggests that EFT is a perfect tool to study the ud 
breakup processes because, once Li^a is fixed by one experiment, predictions can be made 
in other channels such as the solar fusion process pp — > de^u, and supernova short term 
cooling processes such as np —>■ npuu. One experiment which could, in principle, constrain 
Li A is parity-violating ~^d — > enp. The SAMPLE II collaboration |5l| is measuring this 
reaction at MIT-Bates with the intent of studying the strange magnetic moment, iXg. Their 
kinematics are tuned for that purpose, and are not well-suited for an extraction of Li a- It 
would be useful to explore future possibihties that optimize an extraction of the isovector 
axial coupling, Li^a- 

To make a more detailed comparison with potential model calculations, we choose the 
same Li^a values as mentioned above and plot the ratios between different calculations of 
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5 10 15 20 




FIG. 8. Inelastic v{i>)d cross- sections shown as functions of incident v{v) energy. The solid 
curves in the upper graph are KN results while the dot-dashed curves, which lie right on top of the 
solid curves, are NLO results in EFT with Li^a = 6.3 fm^. The solid curves in the lower graph are 
YHH results while the dashed curves, which also lie right on top of the solid curves, are NLO results 
in EFT with Li^a = 1-0 fm^. In both graphs, the highest curves are ud e~^nn cross-sections while 
the middle and lower curves are vd — > vnp and vd — > vnp cross-sections respectively. 

the three i'{v)d reaction channels in Fig. ^ We find that the agreement is excellent between 
our calculations and the potential model results of KN (within 2%), and that we agree with 
the results of YHH to within 10% (the actual agreement is likely better than this - the 
fiuctuations in YHH's calculation may exaggerate the differences somewhat). In making 
these comparisons, we should recall that we have taken into account isospin splitting in 
the ^5*0 rescattering amplitudes for nn and np final states. The details can be found in 
the Appendix. Further, for np final-states, our fit for the ^5*0 channel is not as good as 
for the ^5*1 channel. However, it begins to deviate from the measured phase shifts only 
above relative momenta of 50 MeV. The sensitivity of the calculation to these momenta is 
suppressed by the factor of uj'\k.'\ in the differential cross-section for NC and CC scattering 
(eq. (^OD ). u' oc {p^ax ~ V^) means that the calculation is most sensitive to lower relative 
momenta. 
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FIG. 9. Ratios between different calculations of the inelastic v{i>)d cross- sections shown as 
functions of incident energy. The first graph corresponds to cteft / (^KN while the second and 
third graphs correspond to oeft / c^yhh and oyhh / c^kn respectively. In all graphs, boxes are the 
ud — > i^np ratios, while triangles and diamonds are the vd — > i>np and ud — > e^nn ratios. 

Another interesting quantity is the ratio R between CC and NC scattering cross-sections, 



R 



(83) 



This ratio is zero at threshold {E^ = 4.03 MeV). Above threshold, it is expected to be 
insensitive to the treatment or modeling of short distance physics and, in fact, changes by 
only 5% at NLO above 5 MeV when Li ^ varies over a wide range, from —20 to 40 fm^. 
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FIG. 10. Pd scattering CC to NC cross-section ratios R shown as functions of incident v 
energy. The upper graph corresponds to the NLO EFT result with Li^a = —20 fm^ (boxes) and 
40 fm^ (triangles). The lower graph correspond to the ratios Reft / Rkn (boxes) and Reft/ Ryhh 
(triangles) with input Li^a = 3.7 fm'^ for the evaluation of Reft- 



Choosing Li a = 3.7 fm'^, our values of R agree well with those obtained by KN and YHH 
over the full range of neutrino energies studied (Fig. p!0| ). 

One of the strengths of effective field theory is that it provides us with the ability to 
quantify theoretical uncertainty. It remains for us to quantify the precision of the EFT 
calculations presented here. Naive power counting tells us that the NLO calculation has a 
10% uncertainty from higher order corrections, even if we were able to fit the counterterm 
Li A from experimental data (such as the ~e^d scattering process mentioned earlier). This 
estimate is consistent with the observation that setting Li ^ to +6, 0, and —6 fm'^ (sizes 
suggested from the dimensional analysis of eq. (p2|)) corresponds to the NLO contribution 
being 10, 20 and 30% (respectively) of the LO contribution. As the NLO contribution is 
naively a 30% effect, the NNLO effect could be estimated to be ~ 10%. However, if the sign 
of Li A is positive, as preferred by our comparison to the potential model calculations of 
KN and YHH, the uncertainty from NNLO could be only ~ 4%. This latter scenario would 
also explain why we can fit the potential model results at (in general) the few percent level, 
given that higher order contributions are partially incorporated in very different fashions 



in different calculations. The size of NNLO effects will be studied further ||5^ using EFT 
without pions . 
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VII. CONCLUSIONS 



We have presented analytic differential cross-sections for elastic and inelastic neutrino 
deuteron scattering processes using effective field theory. For elastic scattering, the deuteron 
axial form factor arising from strange matrix elements, and the deuteron strange magnetic 
moment form factor are computed to NLO with two-body current dependence. For inelastic 
scattering, two neutral current processes ud unp, Vd Vnp and one charge current 
process ud — > e~^nn are computed to NLO with an isovector axial two-body matrix element 
whose value is yet to be fixed by experiment. Potential model calculations done by Kubodera 
et al. and Ying et al. are reproduced, with a high degree of accuracy, by choosing different 
values for the two-body matrix element. This implies that the differences between the two 
potential model calculations lie in their treatment of short distance physics. The charged 
current to neutral current Vd cross-section ratio is confirmed to be insensitive to short 
distance physics, and the same ratio is obtained by potential models and this calculation 
within our intrinsic uncertainties (conservatively estimated to be 5%) for the full range of 
incident neutrino energies studied, up to 20 MeV. The two-body matrix element could be 
fixed using the parity-violating process ~e'd — > enp. 

There still remains the need to calculate the other charged current process, vd — > e~pp, 
in EFT. This is the primary reaction channel at SNO and issues of its precision and accuracy 
should not be simply inferred from the processes studied here. The complications introduced 
because of coulomb interactions in the final state make it more appropriate to discuss this 
process elsewhere |Q , along with details of the angular distributions of the charged-current 
reactions 
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APPENDIX 

This appendix summarizes how the A^A^ scattering strong interaction parameters can 
be fit to phase shift data by matching onto the effective range expansion (ERE). The NN 
scattering amplitude A has a power series expansion m. Q; A = A-i + Aq + where the 



24 



subscripts denote the powers in Q. A^i and have been calculated in The standard 
parameterization of the A^A^ scattering amplitude is given by 

A = ^ \ . . (84) 

Mat p cot — ip 

where p is the relative momentum and 6 is the phase shift. The effective range expansion is 
a power series expansion of p cot 6, yielding 

pcotS = ~ + ^^qP^ H , (85) 

where a is the scattering length and tq is the effective range. For the ^5*0 channel, we can 
rewrite eq. ( ^4|) as 



jERE ^ 47r 1 



1 ri'^"V' 

1 J !^ L. L 



_4EH£^^EH£^ , (86) 



where we have performed a Q-expansion. At LO, we can relate the coefficient Cq^° to 
effective range parameters through equating 

A_,{p)=A^Tip) (87) 

yielding 

At NLO, ^0 does not have the same p dependence as A^^'^. The NLO matching is performed 
near p = 0. As p — > 0, 

^r-^(p) ^ -|^a(^^o)V(^^oy _ (89) 

The matching yields 

D^^''\m^) = , (90) 

and 

The ERE parameters 

a„p = -23.714 ± 0.013 fm , ro„^ = 2.73 ± 0.03 fm , 
a„„ = -18.5 ± 0.4 fm , ro„„ = 2.80 ± 0.11 fm , 
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are different between n-p and n-n systems Talcing into account this observed violation 
of charge independence, we find for the n-p system 



^('^"^(mO = -3.56 fm' , D^'^°^(m^) = fm^ , C^"^''^(m^) = 2.79 fm 



while for the n-n system 

d'^°)(m,) = -3.50 fm^ , D^' ^'\m^) = fm^ , d'^^'Vm^) = 2.71 fm^ 



For the ^5*1 system, pcot5 is usually expanded around the deuteron pole, 

pcot5 = -7 + "^^^(p^ + 7^) H , 



(92) 



(93) 



(94) 



such that we can rewrite eq.(|8^ as 

in 1 



A 



ERE 



Mm —7 — ip 

I ERE I A ERE 



after performing a Q-expansion. The LO matching yields 



Mn {fi - 7) 



(95) 



(96) 



At NLO, the matching is performed at p = and 27, so that the scattering length is 
reproduced and the residue of the deuteron pole is not changed in the NLO EFT amplitude. 
That yields two conditions 



TV 



(97) 



Matching at /i = m^, these two conditions can be solved for the coefficients D^^^ and C^^^ 

(98) 



ml {rriT, — 7) 



2 ' 



and 



Numerically we find 



2Df^^)(m.)m^ + ^) + 



9a 



+ 27 
m 1 



(99) 



C^''^^^(m^) = -5.64 fm^ , Z^f ^^^(m^) = 1.46 fm* , C^"''^^(m^) = 10.05 fm' 



(100) 



The fit to the n-p scattering phase shifts is shown in Fig. ^ 
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FIG. 11. n-p scattering phase shifts (measured in degrees) for the "^Si (upper graph) and ^Sq 
channels (lower graph). The dashed curves are the LO results and the solid curves are the NLO 
results from eqs.( 10(\ ,9^). The dotted curves are the results of the Nijmegen partial wave analysis. 
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